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We show that keeping a constant lower limit on the net-time headway is the key mechanism 
behind the dynamics of pedestrian streams. There is a large variety in flow and speed as functions 
of density for empirical data of pedestrian streams, obtained from studies in different countries. The 
net-time headway however, stays approximately constant over all these different data sets. 

By using this fact, we demonstrate how the underlying dynamics of pedestrian crowds, naturally 
follows from local interactions. This means that there is no need to come up with an arbitrary flt 
function (with arbitrary flt parameters) as has traditionally been done. 

Further, by using not only the average density values, but the variance as well, we show how the 
recently reported stop-and-go waves [Helbing et al.. Physical Review E, 75, 046109] emerge when 
local density variations take values exceeding a certain maximum global (average) density, which 
makes pedestrians stop. 

PACS numbers: 89.40.-a, 87.23.Ge, 12.00.00 



I. INTRODUCTION 

With an increasing population and with more cost ef- 
fective transportation, mass gatherings become more fre- 
quent. The total size of such gatherings are often as large 
as millions of people, for example during the inauguration 
ceremony of president Obama |T] and the Hajj pilgrimage 
to Mecca [2]. 

To guarantee the safety of the participants during such 
large mass gatherings, careful planning needs to be car- 
ried out by the organizer. During the last decades, nu- 
merous empirical studies |2][3l|4l|5l|6l[7l|8l|9l[T0l[TT] 
[T2] [13] have been performed on pedestrian crowds in dif- 
ferent countries, in order to understand the dynamics of 
these crowds. Even though an understanding of crowd 
dynamics is a prerequisite for being able to plan a mass 
gathering, there is still no concensus on some of the most 
basic relations, such as how the flow of people (people per 
meter per second) depends on the crowd density (people 
per m^). Misconceptions of these basic relations may 
result in serious safety risks during mass gatherings |13j. 

Let us now start from the bottom up, and show how 
local interactions lead to certain flow-density relation- 
ships for the stream of pedestrians. Since movement and 
avoidance patterns of pedestrians tend to be rather com- 
plex, the traditional way to reduce complexity is to find 
a relation of the flow Q {m~^s~^) as a function of the 
average density g (m~^). Using this relation, called the 
fundamental diagram, has been successful to some extent, 
but unfortunately there are large variations of these re- 
lations, among empirical studies carried out in different 
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countries. All these studies agree on that walking speed 
of pedestrians is a decreasing function of density, but 
they disagree on how this function looks like. We will 
now demonstrate how the net-time headway, as a result 
of finite reaction times, is the key mechanism which can 
explain the discrepancies between data sets from different 
studies. 



II. REACTION TIME 

It is known from traffic science that finite reaction 
times are needed to explain instabilities in traffic fiows 
|14] . For pedestrian-fiow dynamics, the role of finite re- 
action times has not been investigated in detail. By doing 
so, it turns out that the finite reaction time gives rise to 
a certain net-time headway, which is needed as a safety 
headway, to prevent accidental physical encounters with 
surrounding pedestrians. 

Many-particle simulations |151IT6] coupled with empiri- 
cal pedestrian-trajectory data [TJ reveal the probability- 
density function of delay times T'^ from a walking exper- 
iment [l^ where two pedestrian streams are intersecting 
at a 90° angle. The resulting distribution of delay times 
are shown in Fig. [l] 

Interestingly, the probability-density function of the 
delay curve is bi-modal. The first peak occurs at lower 
times than the typical response times, to visual or acous- 
tic cues [19J [23. Therefore, this peak must correspond 
to anticipated movements of the surrounding pedestrians. 
The second peak at around 0.45 s occurs at times which 
are significantly larger than the previously mentioned re- 
sponse times, but also lower than response times involv- 
ing conscious reactions |14| [21]. Therefore, we conclude 
that this second peak corresponds to an unconscious re- 
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Figure 1: Probability-density function of the delay T'*. The 
error bars correspond to one standard deviation. Interest- 
ingly, the probability-density function is bi-modal. When the 
surrounding pedestrians are acting in a way that is easy to 
predict, extrapolation allows to anticipate their behaviors, 
while a delayed reaction results in cases of unexpected be- 
haviors. 



sponse, which is more complex than a simple reaction. 
In fact, it has been shown that reactions where there are 
more than one possible response (choice reaction time) 
as well as reactions to more complex cues (recognition 
reaction time) take significantly longer time |22] . 

We interpret the bi-modality as follows: When the sur- 
rounding pedestrians act in a way that is easy to predict, 
extrapolation allows to anticipate their behaviors, while a 
delayed reaction results in cases of unexpected behaviors. 



III. MODEL 

There has been a rich amount of microscopic models of 
pedestrian dynamics, for example the social-force model 
[151 ITG] and cellular-automata models [23} [231 ISS] ■ These 
models are able to reproduce various self-organization 
phenomena, such as lane and stripe formation [26], freez- 
ing by heating [23, Mexican waves in excitable media 
|28) . intermittent outfiows [29j, stop-and-go waves and 
crowd turbulence |15] . 

When measuring empirical pedestrian fiows and densi- 
ties and then fitting a suitable curve to the data, one ob- 
tains a function which is useful for engineering involved 
in planning of pedestrian facilities. This pragmatic fit 
curve, however, does not provide any insight into the 
mechanisms and dynamics behind the pedestrian inter- 
actions and behaviours, leading to the aggregated data. 

However, when plotting the fundamental diagrams ob- 
tained in various empirical studies (Fig. |3] (top)), one can 
see that each of the curves has a similar parabola-like 



5 




Ql . . . . 1 

0.5 1 1.5 2 

Density (m"^) 

Figure 2: The distance between an arbitrary pedestrian a and 
the closest surrounding pedestrian /3, as a function of global 
(average) crowd density g. The solid line shows the average 
value ± one standard deviation as error bars, and the dashed 
line shows the fit curve The data is from Ref. |18] 



shape. Nevertheless, the curves are quite different from 
one measurement site to another. One question remains 
to be answered: What, if any, is the common underlying 
principle of these curves? 

In an attempt to bridge this knowledge gap, let us come 
back to the issue of reaction times, mentioned before. 
Since pedestrians have a typical reaction time T'^=0.45 
s to unexpected behaviors of surrounding pedestrians, it 
would be natural that they compensate the risk of bump- 
ing into others, by keeping a certain safety time headway 
to the surrounding pedestrians |34] . 

To connect the aggregated density to local interactions, 
let us approximate the mean distance between the center- 
of-masses of a pedestrian a and the closest pedestrian /3 
by d = (dap) ~ where g is the global [35] (average) 

density. Note that this would hold only if the pedestri- 
ans were distributed into a square lattice, but for other 
density distributions, it will serve as a fair approximation 
(see Fig. |2j. 

The net di stance is defined a.s d = d — 2r, where r = 
^/i^\/Bmaxj is the effective radius of a pedestrian, and 
Qmax is the largest measured density. 

Assuming that the predecessor P [36] would sud- 
denly stop it would take T = d/va = {^/y/g — 
^/y/Qmax)/va seconds before a physical encounter with 
pedestrian a occurs, if Va is the speed of pedestrian 
a. We now show how the net-time headway T depends 
on the global density g by applying the above scheme 
to empirical data determined from different studies (see 
Fig.|3];bottom)). 

Note that T saturates at a constant value, that is very 
similar to the response time to unexpected behaviors (see 
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Fig. [ij. However, in the data of Ref. p] there is a transi- 
tion at very high densities, where T suddenly increases. 
This can be interpreted in at least two ways: 

• Hypothesis 1: When the density is very high, pedes- 
trians start to have fear of crushing or asphyxia |30) . 
and therefore want to increase the space around 
themselves (leading to higher net-time headways 
f). 

• Hypothesis 2: If the space in front of a pedestrian is 
too small (or the velocity is too low) it will no longer 
be possible to take normal steps. Rather, pedes- 
trians will completely stop until they have gained 
enough space to make a step. 

In previous work j3^. Hypothesis 1 has been used. In 
this study, however, we will investigate Hypothesis 2. 
This interpretation would naturally explain the empir- 
ically observed stop-and-go waves analyzed in Ref. [2], 
and would further imply: Above an average density of 5 
persons per rn? , the fundamental diagram will no longer 
describe the dynamics of the crowd well, since the flow 
rate is then alternating between movement and standstill 
rather than continuous. 

In an attempt to unify all fundamental diagrams in 
the same framework, the following scheme is proposed: 

Each pedestrian a has a free speed t;° = Vmax (which 
is an upper speed limit, occurring when g ^ Q). Each 
pedestrian also has a lower limit Vmin of the speed. For 
V < Vmin, pedestrians can no longer make normal steps, 
and would rather stop completely. For simplicity, these 
values are assumed to be the same for all pedestrians. 
It has been reported in Ref. |3l] that the free (uncon- 
strained) headways are exponentially distributed, where 
the constrained headways on the other hand, are limited 
by a desired minimum headway. Therefore, we propose 
that the net-time headway T is the key control parameter 
for the fundamental diagram. That is, pedestrians will 
decrease their speeds, if necessary, to assure a constant 
lower limit of the net-time headway T. The fundamental 
diagram can now be specified as: 
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and bounded by [Vmin,Vmaxi 

It has been shown in Ref. |32] that each average density 
g corresponds to a distribution of local densities p, and 
we therefore approximate the density distribution with 
a Gaussian distribution 7V(g, g/3) with mean g and 
standard deviation \/~gJ^ (see Fig. kl. 

According to Hypothesis 2, defineoabove, we get an ex- 
tra constraint, saying that pedestrians will stop walking 
if they are too close to other pedestrians, which happens 
for p > gmax (physical interaction). They will then re- 
sume walking again when they have enough space L for 
taking a step. Since one step (for low walking speed) 
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Figure 3: (Color online) Top: Flow as a function of density, 
for data from a number of empirical studies. Bottom: The 
net-time headways T as a function of density q. T is most 
often bounded by a constant lower value of about 0.5 s. In 
the data of Ref. j2], however, there is a transition for high 
densities where T suddenly increases. The data sets are the 
same as used in Ref. i.e. the data from (Helbing et al.) 
correspond to local densities and flows. 



needs approximately L « 0.5 m [31], we get a new net- 
time headway T' — L/vmin ~ 10 s, whenever p > gmax- 
The fraction of pedestrians that are physically col- 
liding with others, can be measured by integrating the 
probability-density-function of the Gaussian distribution 
(see Fig. |6] (top) ). 



fstop — 



Af{p)dp 



(2) 



with mean p, = g and standard deviation cr = ^/gji. 
Then, the mean net-time headway (see Fig. |6] (bottom)) 
is given by the fraction of stopped pedestrians as 



(T) — (1 — fstop)T + fstop- 



L 



(3) 



Figure [5] shows generated fundamental diagrams from 
Eq. all with the parameters T = 0.5 s, gmax = 5.4 m~^. 
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Figure 4: (Color online) Distributions of local densities p 
for three different global (average) densities q = 1.6,3 and 
5 pedestrians/m'^. The data comes from Ref. For each 
global density g a Beta distribution is fitted (dashed lines). 
However, Gaussian distributions (solid lines) also fit the data 
fairly well. The Gaussian distributions are produced with the 
parameters fi ~ Q and a = ^ q/3. 



and for different values of the free speed -y" = Vmax- Since 
only gives the upper limit of the velocity, fundamen- 
tal diagrams with different converge at high enough 
crowd densities, given that all other parameters are fixed. 
The reason is that, for high density, the movement is 
transformed from individual walking to walking which is 
constrained by other pedestrians. 

We now apply the method outlined above on different 
empirical fundamental diagrams. In all cases we use T — 
0.5 s and Vmin = 0.06 m/s. Starting with Weidmann's 
[3] fundamental diagram, we have the parameters Qmax — 
5.4 m~^ and = 1.34 m/s, which is displayed in Fig. [t] 
(a, b) together with our curve, obtained by Eqs. (fTl to 

Next, we apply our method to the fundamental dia- 
grams of Refs. ID, |2] and [11], and obtain the results 
presented in Fig. [t] (c-h). The fit functions match all 
four different empirical data sets well. All parameters 
are kept constant over the different data sets, except the 
maximum density and the free speed, but these two val- 
ues are obtained from the data, rather than tuned in 
order to fit the data. 
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Figure 5: (Color online) Fundamental diagrams and velocity- 
density relations generated by Eq. ([TJ, for different free speeds 
u". Note how they all converge for large densities. As a 
comparison, the empirical fit curve by Weidmann |3] is shown. 



method is that it follows naturally, without the need of 
an arbitrary fit function. Further, all the parameters are 
measurables, such as the free speed and the maximum 
density. There is not a single free parameter that must 
be tuned in order to fit the different data sets. However, 
it should be mentioned that even though the maximum 
density can be estimated, it can normally not be exactly 
determined from the data. This is addressed in recent 
work |33] that may make it possible to obtain culturally 
dependent parameters, such as the maximum density. 
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Figure 6: Top: The mean net-time headway (T) is obtained 
via the fraction of pedestrians who are physically colliding 
with others, and is therefore stopping and temporarily in- 
creasing their net-time headway. This fraction is obtained 
by integrating over the probability-density-function of the 
local-density distribution, starting at local densities p that 
are higher than the maximum global density q. Bottom: The 
mean net-time-headway (T) (solid line) as a function of the 
global density q. The net-time headway without stopping is 
displayed as a dashed line. 



[1] P. Taylor, How Officials Will Control the Crowds at 
Obama's Inauguration, Popular Mechanics, January 20, 
2009. 

[2] D. Helbing, A. Johansson, and H. Z, Al-Abideen, The 

Dynamics of Crowd Disasters: An Empirical Study, 

Phys. Rev. E 75, 046109 (2007). 
[3] U. Weidmann, Transporttechnik der Fufiganger, ETH- 

Ziirich, Schriftenreihe IVT-Berichte 90 (1993). 
[4] M. Mori H. Tsukaguchi, A new method for evaluation 

of level of service in pedestrian facilities. Transportation 

Research A 21(3), pp. 223-234 (1987). 
[5] J. J. Pruin, Designing for pedestrians: A level-of-service 

concept. Highway Research Record 355, 1-15 (1971). 
[6] K. Ando, H. Ota, and T. Oki, Forecasting the flow of 

people. Railway Research Review 45, 8-14 (1988). (in 

Japanese). 

[7] A. Polus, J. L. Schofer, and A. Ushpiz, Pedestrian flow 
and level of service. Journal of Transportation Engineer- 



ing 109, 46-56 (1983). 
[8] R. A. Smith and J. F. Dickie (eds.). Engineering for 

Crowd Safety. (Elsevier, Amsterdam) (1993). 
[9] K. Still, Crowd Dynamics, PhD Thesis (2000). 

[10] K. Teknomo, Microscopic Pedestrian Flow Characteris- 
tics: Development of an Image Processing Data Collec- 
tion and Simulation Model. PhD Thesis, Japan (2002). 

[11] A. Seyfried, B. Steffen, W. Klingsch, and M. Boltes, The 
fundamental diagram of pedestrian movement revisited, 
J. Stat. Mech. P10002 (2005). 

[12] T. Kretz, A. Griinebohm, and M. Schreckenberg, Exper- 
imental study of pedestrian flow through a bottleneck, J. 
Stat. Mech P10014 (2006). 

[13] A. Schadschneider, W. Klingsch, H. Kluepfel, T. Kretz, 
C. Rogsch, and A. Seyfried, Evacuation Dynamics: Em- 
pirical Results, Modeling and Applications, Encyclope- 
dia of Complexity and System Science, B. Meyers (Ed.) 
(Springer, Berlin, 2008). 



6 



a 1.5 



b 1.5r 




2 4 
Density (m"^ ) 



2 4 
Density (m"' 



Figure 7: (Color online) Fundamental diagrams and velocity- 
density relations generated by Eqs. ([T]) to ([3|, assuming a 
constant net-time headway T = 0.5 s and a minimum ve- 
locity Vmin = 0.06 m/s. Red markers represent empirical 
data and solid lines the theoretically expected relationships, 
a) Flow-density data, and b) velocity-density data by Wei- 
dmann compared to a fundamental diagram generated 
with the parameters Qmax = 5.4 m~^ and Vmax = 1-34 m/s. 
The dashed line shows the result when it is assumed that 
no pedestrians are stopping, i.e. fstop ~ 0. c) Flow-density 
data, and d) velocity-density data by Mori and Tsukaguchi 
Ul, compared to a fundamental diagram generated with the 
parameters Qmax ~ 12 m~^ and Umax = 1-45 m/s. e) Flow- 
density data, and f) velocity-density data from Helbing et 
al. [2], compared to a fundamental diagram generated with 
the parameters Qmax = 9-3 m ^ and Vmax = 0.45 m/s. g) 
Flow-density data, and h) velocity-density data from Seyfried 
et al. jll] , compared to a fundamental diagram generated with 
the parameters Qmax ~ 5.4 m~^ and Vmax ~ 1-35 m/s. 



[14] M. Treiber, A. Resting, and D. Helbing, Delays, inac- 
curacies and anticipation in microscopic traffic models, 
Physica A 360, pp. 71-88 (2006). 

[15] D. Helbing and P. Molnar, Social force model for pedes- 
trian dynamics, Phys, Rev. E 51, pp. 4282-4286 (1995). 

[16] D. Helbing, I. Farkas, and T. Vicsek, Simulating dynam- 



ical features of escape panic. Nature 407, pp. 487-490 
(2000). 

[17] A. Johansson, Data-driven modeling of pedestrian 
crowds, PhD Thesis, TU Dresden (2009). 

[18] S. P. Hoogendoorn, W. Daamen, and P. H. L. Bovy, 
Extracting microscopic pedestrian characteristics from 
video data. Annual Meeting Transportation Res. Board 
Pre-print CD-Rom, (Mira Digital Publishing, Washing- 
ton, D.C) (2003). 

[19] R. S. Woodworth and H. Schlosberg, Experimental Psy- 
chology. Henry Holt, New York (1954). 

[20] A. T. Welford, Choice reaction time: Basic concepts. In 
A. T. Welford (Ed.), Reaction Times. Academic Press, 
New York, pp. 73-128 (1980). 

[21] M. Green, How Long Does It Take to Stop? Methodolog- 
ical analysis of driver perception-brake times. Transport. 
Hum. Factors 2 pp. 195-216, (2000). 

[22] F. C. Donders, On the speed of mental processes. Trans- 
lated by W. G. Koster, 1969. Acta Psychologica, 30, 
pp. 412-431 (1868). 

[23] A. Kirchner, K. Nishinari, and A. Schadschneider, Fric- 
tion effects and clogging in a cellular automaton model 
for pedestrian dynamics, Phys. Rev. E 67, 056122 (2003). 

[24] C. Burstedde, K. Klauck, A. Schadschneider, and J. Zit- 
tartz. Simulation of pedestrian dynamics using a two- 
dimensional cellular automaton, Physica A, 295(3-4), 
pp. 507-525 (2001). 

[25] W. G. Weng, T. Chen, H. Y. Yuan, and W. C. Fan, 
Cellular automaton simulation of pedestrian counter flow 
with different walk velocities, Phys. Rev. E 74, 036102 
(2006) 

[26] D. Helbing, L. Buzna, A. Johansson, T. Werner, Self- 
organized pedestrian crowd dynamics: Experiments, sim- 
ulations, and design solutions, Transp. Sci. 39, pp. 1-24 
(2005). 

[27] D. Helbing, I. Farkas, and T. Vicsek, Freezing by heat- 
ing in a driven mesoscopic system, Phys. Rev. Lett, 84, 
pp. 1240-1243 (2000). 

[28] I. Farkas, D. Helbing, and T. Vicsek, Mexican waves in 
an excitable medium. Nature 419, pp. 131-132 (2002). 

[29] D. Helbing, A. Johansson, J. Mathiesen, H. M. Jensen, 
and A. Hansen, Analytical approach to continuous and 
intermittent bottleneck flows, Phys. Rev. Lett. 97, 
168001 (2006). 

[30] W. Yu and A. Johansson, Modeling crowd turbulence 
by many-particle simulations, Phys. Rev. E 76, 046105 
(2007). 

[31] S. P. Hoogendoorn, W. Daamen, Pedestrian behavior at 
bottlenecks, Transp. Sci. 39(2), pp. 147-159 (2005). 

[32] A. Johansson, D. Helbing, H. Z. A-Abideen, S. Al-Bosta, 
From crowd dynamics to crowd safety: A video-based 
analysis. Advances in Complex Systems 11(4), pp. 497- 
527 (2008). 

[33] U. Chattaraj, A. Seyfried, and P. Chakroborty, Compar- 
ison of pedestrian fundamental diagram across cultures. 
Advances in Complex Systems 12, 393 (2009). 

[34] Note that we use the peak at longer delays since this peak 
corresponds to unexpected events. Even though there is 
a peak at lower delay times, there is no guarantee that 
the behaviours of others can always be anticipated, and 
therefore the maximum peak must be used as a safety 
time. 

[35] The global density q is deflned as the number of people 
within a certain area, divided by that area. The local den- 



sity p [2] on the other hand, is defined via a bell-shaped 
weight function, where the influence of close pedestrians 
is larger than the influence of remote pedestrians. 
[36] This assumes a following behaviour for uni-directional 
streams, which is most often the case for empirical pedes- 
trian studies. 



